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Abstract 

We present an implementation of a parton shower algorithm for hadron colliders and electron- 
positron colliders based on the dipole factorisation formulae. The algorithm treats initial-state 
partons on equal footing with final-state partons. We implemented the algorithm for massless 
and massive partons. 



1 Introduction 



Event-generators like Pythia [1,2], Herwig [3,4] or Sherpa [5] are a standard tool in high energy 
particle physics. In these tools the physics of particle collisions is modelled by a simulation with 
different stages - hard scattering, parton showering, hadronisation - to name the most important 
ones. The hard scattering process is calculable in perturbation theory. The same holds - in 
theory at least - for the parton showering process, the relevant scales are still large enough for 
perturbation theory to be applicable. In practise however, one is forced into approximations due 
to the large parton multiplicities. These approximations are derived from the behaviour of the 
matrix elements in singular regions. The matrix elements become singular in phase space regions 
corresponding to the emission of coUinear or soft particles. The first showering algorithms started 
from the collinear factorisation of the matrix elements and approximated colour interference 
effects through angular ordering [6,7] . An exception is the algorithm implemented in Ariadne [8- 
12], which is based on a dipole cascade picture. Most shower algorithms are in the collinear limit 
accurate to the leading-logarithmic approximation. Extensions to the next-to-leading logarithmic 
approximation have been studied in [13-16]. 

Recent years have witnessed significant developments related to shower algorithms, includ- 
ing procedures to match parton showers to fixed-order tree-level matrix elements [17-21] and 
methods to combine parton showers with next-to-leading order matrix elements [22^5]. The 
shower algorithms in Pythia, Herwig and Ariadne have been improved [46^8] and new pro- 
grams like the shower module Apacic++ [49, 50] of Sherpa have become available. Other im- 
provements include the study of uncertainties in parton showers [51-53], as well as showers in 
the context of the soft-coUinear effective theory [54] . 

Of particular importance is the matching of parton showers with next-to-leading order ma- 
trix elements. The pioneering project MC@NLO [31,55-58] used an existing shower program 
(Herwig) and adapted the NLO calculation to the shower algorithm, at the expense of sacrificing 
the correctness in certain soft limits. It is clear that a better but more labour-intensive approach 
would adapt the shower algorithm to NLO calculations. Nowadays in NLO computations the 
dipole subtraction method [59-63] is widely used. Nagy and Soper [35,36] proposed to build a 
shower algorithm from the dipole subtraction terms. 

In this paper we report on an implementation of a shower algorithm based on the dipole 
formalism as suggested by Nagy and Soper. We take the dipole splitting functions as the splitting 
functions which generate the parton shower. In the dipole formalism, a dipole consists of an 
emitter- spectator pair, which emits a third particle, soft or collinear to the emitter. The formalism 
treats initial- and final-state partons on the same footing. In contrast to other shower algorithms, 
no distinction is made between final- and initial-state showers. The only difference between 
initial- and final-state particles occurs in the kinematics. In the implementation we have the four 
cases final-final, final-initial, initial-final and initial-initial corresponding to the possibilities of 
the particles of the emitter-spectator-pair to be in the initial- or final-state. All four cases are 
included, therefore the shower can be used for hadron colliders and electron-positron colliders. 
We implemented the shower for massless and massive partons. Initial-state partons are however 
always assumed to be massless. We use spin-averaged dipole splitting functions. The shower 
algorithm is correct in the leading-colour approximation. As the evolution variable we use the 
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transverse momentum in the massless case, and a variable suggested in [47, 64] for the massive 
case. The variable for the massive case reduces to the transverse momentum in the massless 
limit. Schumann and Krauss report on a similar but separate implementation of a parton shower 
algorithm based on the dipole formalism [65]. 

This paper is organised as follows: In section 2 we review basic facts about the colour de- 
composition of QCD amplitudes and the dipole formalism. In section 3 we discuss the shower 
algorithm. In section 4 we present numerical results from the parton shower simulation program. 
Finally, section 5 contains the summary. Technical details can be found in the appendix. Ap- 
pendix A discusses the case of a massless final-state emitter and a massless final-state spectator 
in detail. Appendix B describes the construction of the four-momenta of the {n+ 1) -particle state 
in all cases. This appendix is also useful in the context of a phase-space generator for the real 
emission part of NLO computations. 

2 QCD amplitudes and the dipole formalism 

In this section we briefly review the colour decomposition of QCD amplitudes and the dipole 
formalism. 

2.1 Colour decomposition 

In this paper we use the normalisation 

Trr^r^ = l5«* (1) 

for the colour matrices. Amplitudes in QCD may be decomposed into group-theoretical factors 
(carrying the colour structures) multiplied by kinematic functions called partial amplitudes [66- 
70]. The partial amplitudes are gauge-invariant objects. In the pure gluonic case tree level 
amplitudes with n external gluons may be written in the form 



Jln{l,2,...,n) = 




where the sum is over all non-cyclic permutations of the external gluon legs. The quantities 
A„ ( o 1 , . . . , o„) , called the partial amplitudes, contain the kinematic information. They are colour- 
ordered, e.g. only diagrams with a particular cyclic ordering of the gluons contribute. The choice 
of the basis for the colour structures is not unique, and several proposals for bases can be found 
in the literature [71,72]. Here we use the "colour-flow decomposition" [72,73]. This basis 
is obtained by replacing every contraction over an index in the adjoint representation by two 
contractions over indices / and j in the fundamental representation: 

ya^a ^ y«5«^^^ = y« {^T-jTj'^ = (|V27;.^V«) (^27,7^^) . (3) 
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As a further example we give the colour decomposition for a tree amplitude with a pair of quarks: 

^„+2{q,l,2,...,n,q) = (^-^ Y,^i^M^i^ijo2--^io„j9^n+2{q,Ci,C2, ...,Cn,q), (4) 

where the sum is over all permutations of the gluon legs. The tree amplitude with a pair of 
quarks, n gluons and an additional lepton pair has the same colour structure as in eq. (4). In 
squaring these amplitudes a colour projector 

has to applied to each gluon. In these examples we have two basic colour structures, a colour 
cluster described by the "closed string" 

and a colour cluster corresponding to the "open string" 

Bom amplitudes with additional pairs of quarks have a decomposition in colour factors, which 
are products of the two basic colour clusters above. The colour factors in eq. (2) and eq. (4) are 
orthogonal to leading order in l/Nc- 



2.2 The dipole formalism 

The starting point for the calculation of an observable O in hadron-hadron collisions in perturba- 
tion theory is the following formula: 

{O) = /<i.2/(«)^(^pj^;^ 

j d(?n{Pl,P2;P3,---,Pn+2)0{pi,...,Pn+2)\^n+2\^ ■ (8) 

This equation gives the contribution from the n-parton final state. The two incoming particles are 
labelled pi and p2, while pj, to Pn+2 denote the final state particles. f{x) gives the probability 
of finding a parton a with momentum fraction x inside the parent hadron h. A sum over all 
possible partons a is understood implicitly. 2K{s) is the flux factor, 1 / (27i + 1) and 1 / (272 + 1) 
correspond to an averaging over the initial helicities and ni and n2 are the number of colour 
degrees of the initial state particles, d^n is the phase space measure for n final state particles, 
including (if appropriate) the identical particle factors. The matrix element |^n+2p is calculated 
perturbatively. At leading and next- to-leading order one has the following contributions: 

{0)^0 ^ Jonda^, 

n 

{OfLO ^ J On+ida^ + JOndo^ + JOndo^. (9) 

n+1 n n 
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Here we used a rather condensed notation, do^ denotes the Bom contribution, while da^ de- 
notes the real emission contribution, whose matrix element is given by the square of the Bom 

amplitudes with (« + 3) partons dc^ gives the virtual contribution, whose matrix el- 

ement is given by the interference term of the one-loop amplitude ^li^^2 ^^^^ (" + ^) P^o^s 
with the corresponding Born amplitude ^j^2- denotes a coUinear subtraction term, which 
subtracts the initial-state coUinear singularities. Within the subtraction method one constructs an 
approximation term with the same singularity stmcture as do^. The NLO contribution is 
rewritten as 



{Of^^ = J (^On+ida^-Ondc^^+ 1 (Ondo^ + Ondo^ + Onda^^, 



(10) 



n+1 



such that the terms inside the two brackets are separately finite. The matrix element correspond- 
ing to the approximation term dc^ is given as a sum over dipoles [59-63]: 



pairs i,j k^ij 



pairs i J j k^j j 



(11) 



In eq. (11) the labels /, j and k denote final-state particles, while a and b denote initial-state 
particles. The first term describes dipoles where both the emitter and the spectator are in the 
final-state. Dfj denotes a dipole where the emitter is in the final-state, while the spectator is in 

the initial-state. The reverse situation is denoted by 2)^'^: Here the emitter is in the initial-state 
and the spectator is in the final-state. Finally, 'D^J''^ denotes a dipole where both the emitter and 
the spectator are in the initial-state. The full complexity is only needed for hadron colliders; 
for electron-positron annihilation the subtraction terms inside the square bracket are absent. The 
dipole subtraction terms for a final-state emitter-spectator pair have the following form: 

^ ^(0)*/ ~ ~ \{—^k-T^ij) yij,k ^(0) / - _ \ /nx 

(Dij,k = K+2 {ph...,Piij),-,Ph-) 2 9„ ^ K+2{Ph-,Piij),-,Ph-)-a2) 

The structure of the dipole subtraction terms with initial-state partons is similar. Here T, denotes 
the colour charge operator for parton / and Vij^k is a matrix in the spin space of the emitter 
parton (ij). In general, the operators lead to colour correlations, while the Vij^s lead to spin 
correlations. The colour charge operators T, for a quark, gluon and antiquark in the final state 
are 

quark: R* {...qi...) {Tf^) {...qj...) , 

gluon: J?*(.../...)(//-^)^ (.../. 

antiquark : {...qi...) {-Tf^ Jl {-^j-) ■ (13) 

The corresponding colour charge operators for a quark, gluon and antiquark in the initial state 
are 

quark: Jl* {...qi...) {-Tjl) R {...qj...) , 
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gluon: ^* (.../...) (/r^)^ (.../. 
antiquark : ^* (...<?/...) {T^f) A {...qj...) . (14) 

In the amplitude an incoming quark is denoted as an outgoing antiquark and vice versa. 

In this paper we neglect spin-correlations and work to leading-order in l/Nc- Therefore we 
replace the splitting functions Vij^^ by the spin-averaged splitting functions: 

Vij,k ^ {Vij,k) (15) 

In the leading-colour approximation we only have to take into account emitter-spectator pairs, 
which are adjacent inside a colour cluster. For those pairs we obtain for the colour charge opera- 
tors 

(-Tfc-T,-y) ^ f 1/2 emitter (0") is a gluon, 

T?. \ 1 emitter {ij) is a quark or antiquark. 

We introduce the notation 



The functions !P will govern the emission of additional particles in the shower algorithm. The 
spin-averaged dipole splitting functions {V) can be found in [59, 63]. The Heavyside theta- 
functions ensure that the functions will be non-negative. They are needed for splittings be- 
tween an initial- and a final-state particle, since the dipole splitting functions (Vfj) and {V^'') 
may take negative values in certain regions of phase space. In addition, the spin-averaged dipole 
splitting functions for massive partons are slightly modified: Terms related to the soft singularity 
are re-arranged between the two dipoles forming an antenna, in order to ensure positivity of the 
individual dipole splitting functions in the singular limit. 



3 The shower algorithm 

In this section we describe the shower algorithm. We first discuss the colour treatment in sec- 
tion 3.1. The shower algorithm for massless final-state partons is discussed in section 3.2. The 
necessary modifications for initial- state partons are discussed in section 3.3. Finally, massive 
partons are discussed in section 3.4. 

3.1 Colour treatment 

Before starting the parton showers, the partons from the hard matrix element have to be assigned 
to colour clusters. For the simplest matrix elements, like e~^e~ — > qq, the choice is unique: The 
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quark-antiquark pair forms a colour cluster. For the parton shower we work in the leading-colour 
approximation. In the leading-colour approximation we have to take into account only emitter- 
spectator pairs, which are adjacent inside a colour cluster. We have implemented two options: 
In the first one, which we call the "strict leading-colour approximation", we take exactly the 
terms which are leading in an expansion in 1 /Nc and only those. As a consequence, all splittings 
g ^ qqdXQ ignored, as they are colour-suppressed compared to g — > gg. In this approximation 
Cp is replaced by 

Cf ^ \. (18) 

For the second option, which we call the "modified leading-colour approximation", we include 
the splitting g ^ qq and keep Cp as (A/^ — \ )/l/Nc. In this case, if a gluon in a closed string 
splits into a quark-antiquark pair, the closed string becomes an open string. If a gluon in an open 
string splits into a quark-antiquark pair, the open string splits into two open strings. 



3.2 The shower algorithm for massless final-state partons 

We first describe the shower algorithm for electron-positron annihilation. The extension to 
initial-state partons is treated in section 3.3. For the shower algorithm we use as an evolution 
variable 

t = In— ^, (19) 

Q 

where is a fixed reference scale and kj_ is the transverse momentum of a splitting. During the 
shower evolution we move towards smaller (more negative) values of t. We start from a given 
n-parton configuration. In the dipole formalism, emission of additional partons is described by 
an emitter-spectator pair. In the leading colour approximation emitter and spectator are always 
adjacent in the cyclic order. The probability to evolve from ti to t2 (with ti > t2) without any 
resolvable branching is given by the Sudakov factor. For the algorithm considered here, the 
Sudakov factor is given as a product of factors corresponding to the no-emission probabilities for 
individual dipoles' emissions: 

HhM) = YlAikithti). (20) 

i,k 

If parton i can emit different partons. A- j(?i,?2) factorises in turn into different contributions: 

hkithti) = llAij,k{tut2), (21) 

j 

An example is the possibility of a gluon to split either into two gluons or into a ^^-pair. We 
denote the emitter before the splitting by i, while the emitter after a splitting is denoted by /. This 
notation takes into account that the emitter might change its "flavour" due to a splitting, like in 
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the case ofag^qq splitting. Aij^k{ti,t2) is the probabiUty that the dipole formed by the emitter 
/ and spectator k does not emit a parton j. It is given by 



^ij,k{tht2) = exp - j dtCjl j d^unres^ (t - Tjlj 



ii,k 1 



(22) 



where Cji is a colour factor. In the leading colour approximation this factor is non-zero only if i 
and k are adjacent in a colour cluster. Then Cji is obtained from eq. (16) and given by 



^Ik — 



\ for i = g, 
1 for i = q, q. 



(23) 



The dipole phase space is given by 



167l2 



1 z+M 
Jk / dz 



1 



z-(k) 



1- 



K 



4z(l-z)V 4z(l-z) 



(24) 



with 



z±(k) = i(i±v^r^). 

The variable K is proportional to the transverse momentum of the splitting 



(25) 



K = 4 



{Pi + Pk) 



-^2• 



(26) 



T- ^ depends on the dipole invariant mass (/?,■+ pif' and the phase space variable K for the emis- 
sion of an additional particle and is given by 



h,k - q2 



(27) 



With the help of the delta-function we may perform the integration over K, while keeping the 
integration over t and z. Then 



2j 



Kit) 



(Pi+Pk) 



~\2- 



(28) 



fPij^k is the dipole splitting function. As an example we quote the splitting function for the q^qg 
splitting: 



Cf 



{Pi+PkYy 



i-z(i-y) 
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-(1+^) 



K(0 



4z(l-z) 



(29) 



a, is evaluated at the scale = —k\ = ^{pj+ Pif'- The probability that a branching occurs at 
t2 is given by 

L L ^Ik j d^unres^ (?2 - Tjlj ^ij^k^h , ^2) • (30) 
i,k J 

We can now state the shower algorithm. Starting from an initial evolution scale ti we proceed as 
follows: 

1. Select the next dipole to branch and the scale t2 at which this occurs. This is done as 
follows: For each dipole we generate the scale t2,ij,k of the next splitting for this dipole 
from a uniformly distributed number n^ij^^ in [0, 1] by solving (numerically) the equation 

^ij,k{tht2,ij,k) = n,ij,k- (31) 

We then set 

t2 = max (^2,(7,*:) . (32) 

The dipole which has the maximal value of ?2,j7,A: is the one which radiates off an additional 
particle. 

2. If t2 is smaller than a cut-off scale tmin, the shower algorithm terminates. 

3. Next we have to generate the value of z. Again, using a uniformly distributed random 
number r2 in [0, 1] we solve 

Z Z+{t2) 

j dz'j{t2,z!)'Pij,k = ^2 j dz!J{t2,z!)'Pij,k, (33) 

Z-{t2) Z-{t2) 

where the Jacobian factor J{t2,z) is given by 

4. Select the azimuthal angle {|). Finally we generate the azimuthal angle from a uniformly 
distributed number in [0, 1] as follows: 

^ = 27Cr3. (35) 

5. With the three kinematical variables t2, z and (|) and the information, that parton / emits a 
parton 7, with parton k being the spectator, we insert the new parton j. The momenta pj 
and pi of the emitter and the spectator are replaced by new momenta pi and pk- The details 
how the new momenta pi, pj and pk are constructed are given in the appendix B. 
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6. Set ti = t2 and go to step 1. 

Remark: Step 1 of the algorithm is equivalent to first generating the point t2 from a uniformly 
distributed number ri in [0, 1] by solving (numerically) the equation for the full Sudakov factor 



and then selecting an individual dipole with emitter i, emitted particle j and spectator k with 
probability [74] 



3.3 The shower algorithm with initial-state partons 

In this subsection we discuss the necessary modifications for the inclusion of initial-state partons. 

In the presence of initial-state partons there is no separation into final-state showers and initial- 
state showers. Initial-state radiation is treated on the same footing as final-state radiation. The 
algorithm generates initial-state radiation through backward evolution, starting from a hard scale 
and moving towards softer scales. Therefore the shower evolves in all cases from a hard scale 
towards lower scales. 

Final-state emitter and initial-state spectator 

For an initial-state spectator we modify the Sudakov factor in eq. (22) to 



where Xa is the momentum fraction of the initial hadron carried by a, while Xa is the momentum 
fraction carried by a. The initial parton of the n-particle state is denoted by d, while the initial 
parton of the (n -|- 1) -particle state is denoted by a. We set 



A(?i,?2) = n, 



(36) 




(37) 




(38) 



X = 



(39) 



The unresolved phase space is given by 




(40) 



The transverse momentum between / and j is expressed as 



-kl = ^-^z{l-z){-2pjpa) 



(41) 
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and Tjg^ is therefore given by 

A subtlety occurs for the emission between a final-state spectator and an initial-state emitter. We 
discuss this for the splitting q — > qg. The spin-averaged splitting function for the q^qg splitting 
is given by 

(y^) = 87ta,C^[^-^-^^-(l+z)j. (43) 

In contrast to the final-final case this function is not a positive function on the complete phase- 
space. It can take negative values in certain (non-singular) regions of phase-space. This is no 
problem for its use as a subtraction terms in NLO calculations, but prohibits a straightforward 
interpretation as a splitting probability for a shower algorithm. However, since negative val- 
ues occur only in non-singular regions, we can ensure positiveness by modifying the splitting 
functions through non-singular terms. The simplest choice is to set 

%. = ji^'l^im- (44, 

For a final-state emitter we eliminate the jc- integration with the help of the delta- function: 
1 



Xa ^+ K(f) ^ + 4z(l-z) ^ ^'^^^ 



(45) 



For the boundaries we obtain 

K(0<^, z-{t)<z<z+{t), z±(0 = ^l±^l-K(0^). (46) 

The modifications to the shower algorithm are as follows: The dipoles for the emission from a 
final-state emitter with an initial-state spectator are included in the Sudakov factor in eq. (20). 
With this modification steps 1 and 2 are as above. Let us define 

(1, if / and n are final- state particles, 

x"/(x°'f) ' if / = <5[ is an initial-state particle, (47) 



xbf[xb,t) ^ if /J = ^ is an initial-state particle and / is a final-state particle. 



In step 2 we replace formula (33) by 

j dz'J{t2,z')fij,a^ij,a = j dz'J{t2,z')fij,a^ij,a, 

Z- 

with the Jacobian 



Z+(t2) 

I 

Z-ih) Z-{t2) 



(48) 



1 + 



Steps 4 to 6 proceed as in the case described above. 
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Initial-state emitter and final-state spectator 

For an initial-state emitter a. with a final-state spectator / the Sudakov factor is given by 



(50) 



r.r = in-^^ — ^ ' — — (5T) 

^a,i • ^-^"^^ 



The unresolved phase space is again given by eq. (40). The transverse momentum between a and 
j is given by 

-kl = ^-^il-z){-2pjpa) (51) 

and Tqj is given by 

{-2pjpa)_il-x){l-z) 

xQ 

For a initial-state emitter we eliminate the z-integration with the help of the delta-function: 

/■ , _ / _ \ Kit) X , K(t) X , , AQ^e^ 

y<(z8(r-7i,)=Ai_ = (53, 

For the boundaries we obtain 

1 —x~ 1 

K(f)<4 JC<JC+(0, x^{t) = -r^. (54) 



X 



a 1 + 



K(0- 

4 



There are no new modifications to the shower algorithms compared to the case for a final-state 
emitter and an initial-state spectator, except that in step 3 we now generate the value of x accord- 
ing to 

y" dx!j{t2,x!)faj^i'Pai,i = ^2 J dx'j{t2,x')faj,i'Paj,i, (55) 

Xa Xa 

with the Jacobian 

= (56) 

Initial-state emitter and initial-state spectator 

For an initial-state emitter a with an initial-state spectator b the Sudakov factor is given by 

/ 'l 



Kj,b{tht2) = exp 



/*Q,/.w.5(,-r,,)|^.„,,). (57, 

\ '2 
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In this case we do not rescale the momentum of the spectator, but transform all final-state mo- 
menta. Therefore no factor 



Xbf{xb,t) 

appears in the Sudakov factor. The unresolved phase space is given by 

1 l-x 



Xa 

The transverse momentum between a and j is given by 



X 



and T~ ^ is given by 



{^PaPh) (1 -x)v 



We integrate over v with the help of the delta-function: 

\-x 



For the boundaries we obtain 



with the Jacobian 



(58) 



-kl = ^ '-v{2pap~^) (60) 



T^fi = ^""^^^^i (61) 



/■ . o / \ K(n X K(t) X , , 42V 

/.v8(,-7V,)=Ai_ .= Ai_K(,)^^. (62) 



K(»)<4^i-^, x<x+(f), x+W = ^2 + !^-^K«)+^|^l. (63) 
In step 3 of the shower algorithm we again select x according to 

X X+(f2) 

^ dx'j{t2,x')faj^b'Paj,b = f dx'j{t2,x')faj^b^aj,b, (64) 



= (65) 
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3.4 The shower algorithm for massive partons 

In this subsection we discuss the modifications of the shower algorithms due to the presence of 
massive partons. We first address the issue of a splitting of a gluon into a heavy quark pair. This 
mainly concerns the splitting of a gluon into fe-quarks. We will always require that initial- state 
particles are massless. Therefore for processes with initial-state hadrons we do not consider 
g QQ splittings. Calculations for initial-state hadrons should be done in the approximation of 
a massless ^-quark. In the case of electron-positron annihilation the parton shower affects only 
the final state. Here we can consistently allow splittings of a gluon into a pair of massive quarks. 
As evolution variable we use in the massive case 

t = \n ^ ' ' (66) 

Q 

This choice reduces to eq. (19) in the massless limit and is suggested by dispersion relations for 
the running coupling [47,64]. 

Final-state emitter and final-state spectator 

The unresolved phase space is given by 

\2 



% = n2 • (68) 



fdi^unres = ^-^^^^ {l - - - df [Hh ^k)]' ' J dj {I - j) J dz, (67) 

where the reduced masses im and the boundaries on the integrations are defined in appendix B in 
eqs. (1 19)-(122). T- -^ is given by 

{{Pi + Pkf -inj-mj- ml^ yz{ 1 - z) 

a 

Again, we have to ensure that the splitting functions are positive. The original spin-averaged 
dipole splitting functions can take negative values in certain regions of phase-space. In the mas- 
sive case the negative region can extend into the singular region. The problem is related to the 
soft behaviour of the dipole splitting functions. Since a squared Bom matrix element is positive 
in the soft gluon limit, the negative contribution from a particular dipole is compensated by the 
contribution from the dipole, where emitter and spectator are exchanged. The sum of the two 
contributions is positive in the singular region. Therefore we can cut out the negative region from 
the first dipole and add it to the second dipole. The second dipole will stay positive. 
As in the massless case we eliminate the y-integration: 

y+ Z+iy) Zmax 

fdyil-y) I dzb(t-Tj^ = fdzyil-y), (69) 

y- Z- (y) Zmin 



4z(l— z)' {pj+ pj^)^ — mf — mj —m 



T 
k 



14 



The physical region is defined by 



i'-W^i) [h^'-'^'-^'-'H-i-W^] >- 0. (70) 



with 



% = for/eO,M}. (71) 

This equation is solved numerically for Zmin and Zmax- Then z is generated according to 

Z Zmaxitl) 

J dz'J{t2,z')^ij,k = ''2 f dz'J{t2,z')^ij,k, (72) 

Zminitl) Zminih) 

with the Jacobian 

7(r,z) = (73) 

Final-state emitter and initial-state spectator 

The unresolved phase space is given by 

XS Z-{x) Z-iXa) Xa 

The integration boundary is given by 

\-x{\-n^) /u^ + z{l-n^) \2piPa\ 

Tin is given by 



_ -/^i + (1 _ (-2py,) (1 -^)z(l -z) 



For a final-state emitter we eliminate the jc- integration with the help of the delta- function 
1 



Kit) ^ + 4z(l-z) ^ "^'^"^ 



^ (77) 
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For the boundaries we obtain 



The boundary on x{t) is given for /i^ < (1 —xg) /xq by 



Xn 



K{t) < 



1 -Xa 



For (l—Xa) jxa < we have 



k(0 < 



1 -Xa 



Xn 



1- 



< 


Xa 


1- 


^\ 


1 + 


^ / 



z is generated according to 

Z Z+{t2) 

j dz'J{t2,z')fij,a'Pij,a = f dz'J{t2,z')fij,a^> 

Z-{t2) 



i],ai 



with the Jacobian 



1 



1 + 



4z(l-z) 

K(f) 



Initial-state emitter and final-state spectator 

% is given by 

{-2p.pa){\-x){\-z) 



Taj = In- 



xQ^ 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



For an initial- state emitter we eliminate the z-integration with the help of the delta-function: 



1 

/ dzd{t-T,^.) 



For the boundaries we obtain 



<t) X „ _ , k(0 X 4(2V 

4 (1-x)' 4 (1-x)' ^^'^ {-IpjpaY 



(84) 



K(0< 



4(1 



Xa[^ -Xa{l -fl^ 



X<X4-(?), X+(?) 



2 + - 
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+//2K(r) 



2(l + ]f (1-^2)) ' 



(85) 
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The value of x is generated according to 



X+{l2) 



J dx'j{t2,x')faj,i^aj,i ^ r2 J dx'j{t2,x)faj,i^aj,i, (86) 



with the Jacobian 



= (87) 



4 Numerical results 



In this section we show numerical results obtained from the parton shower. We first discuss 
in section 4.1 observables related to electron-positron annihilation and then in section 4.2 the 
shower in hadron collisions. The shower algorithm depends on two parameters, the strong cou- 
pling and the scale Qmin- For the strong coupling we use the leading-order formula 

An 2 
a,(/i) = ^Q=n--Nf. (88) 

The cut-off scale Qmin gives the scale at which the shower terminates. As our shower is correct 
in the leading-colour approximation, we also study the effects of different treatments of sub- 
leading colour contributions. As described in section 3.1 we have implemented two options: The 
strict leading-colour approximation and the modified leading-colour approximation. Numerical 
differences from these two options will give an estimate of uncertainties due to subleading-colour 
effects. 



4.1 Electron-positron annihilation 

For electron-positron annihilation we use (Xs{mz) = 0.118 corresponding to A5 = 88 MeV. We 
start the shower from the 2 — > 2 hard matrix element e'^e" — > qq. We first study the event shape 
variables thrust, the C-parameter and the D-parameter. The distributions of the first moments 
of these observables are shown in figure 1 for two choices of the cut-off parameter: Qmin = 
1 GeV and Qmin — 2 GeV. The distributions are normalised to unity. The different prescriptions 
for the colour-treatment do not change the distributions significantly. In figure 2 we show the 
distributions for the four-jet angles. Again we start from the 2 — > 2 hard matrix element. The 
particles in an event are first clustered into jets, defined according to the Durham algorithm [75] 
with ycut = 0.008 and the ^-scheme for the recombination. Then events with exactly four jets 
are selected. We consider the modified Nachtmann-Reiter angle [76], the Komer-Schierholz- 
Willrodt angle [77], the Bengtsson-Zerwas angle [78] and the angle between the jets with 
the smallest energy [79]. In the plots we show the results from the different options for the colour 
treatment for Qmin = 1 GeV. A variation of the cut-off scale does not change the distributions 
significantly. 
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4.2 Hadron colliders 



For the Tevatron and the LHC we study Z/y* -production. We start from the 2 — > 2 hard matrix 
element qq — > Z/y* — > /+/~. As parton distribution functions we use the CTEQ 6L1 set [80, 81]. 
For consistency we use here as (mz) = 0.130 corresponding to A5 = 165 MeV. The centre-of- 
mass energy we set to ^/s = 1.96 TeV for the Tevatron and to = 14 TeV for the LHC. We 
require a cut on the invariant mass of the lepton pair of 

mi+i- > 80GeV. (89) 

As cut-off parameter for the parton shower we use Qmin = 1 GeV. In figure 3 we show the 
transverse momentum distribution and the rapidity distribution of the lepton pair for the Tevatron 
and the LHC. 



5 Summary 

In this paper we presented an implementation of a shower algorithm based on the dipole for- 
malism. The formalism treats initial- and final-state partons on the same footing. The shower 
can be used for hadron colliders and electron-positron colliders. We also included in the shower 
algorithm massive partons in the final state. We studied numerical results for electron-positron 
annihilation, the Tevatron and the LHC. 
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A Sudakov factors for massless final-state partons 

In this appendix we discuss in more detail the Sudakov factors for massless final-state partons. 
This case is simple enough that one integration can be done analytically. The spin-averaged 
dipole subtraction terms in four dimensions are 



with 



To 



g^qq 



= Cf 
= Ca 



8:ra.v (/./-) 1 
Sijk y 

^nasifP') 1 
Sijk y 



l-z{l-y) 
2 



+ 



l-zil-y) l-(l-z)(l-3;) 



-4 + 2z(l-z) 



= r.^^i[i-2z(i- 
Sijk y 



(90) 



Sijk = {pi+Pj+Pkf = {Pi+PkY 



(91) 
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The dipole phase space measure is 



1 z+(k) 
fd. f dz- ' 



"""^ " l6n^J"-^J ^^4zil-z)V 4z(l-z) 

z_(k) 



1- 



K 



with 



The strong coupUng is evaluated at the scale fi^ = —k']_ : 
The Sudakov factor is given by 

^ij,kitl , h) = exp - dtCi ~k J di^unres^ {t - Tjij fP, 
\ h 

For the splitting q ^ qgwe obtain 



ij,k 



K+ Z+(K) 



A,,,,(ri,r2) = exp -C^C^/^^ / ^z(l-,) 



Z-(K) 



K_=4— e'^ K+ = min( 1,4— e^i ) , y = 



1-zil-y) 



with 

^ijk \ ^ijk 

The integration over z can be done analytically 

2 



K 



4^(1^' ^ " 4^'^'^- 



fdzil-y) 



= -^z^-z+|[lnz-21n(l-z) 



4 + K 



.1-^(1 -y) 

^Klnz + ln (k + 4(1 -z)2) + \/Karctan ^-^(1 -z)^ 



The same holds for the other splittings. Therefore we obtain for the Sudakov factors 



A / N ; ^ fdKas{lKSijk) 



2ti 



{l'ij,ki^,Z+)-'Pij,ki^,Z-)) 
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where C is a colour factor and equal to 

C = 

The functions I'ij^k (k, z) are given by 



Cp for q^qg, 
Ca for g^gg, 
Tr for g^qq. 



-iz2_^+|[ln^_21n(l-z)] 



4 + K 



^Klnz + ln (k + 4(1 -zf) - \/Karctan ^-^(1 -z)^ 



232.1 ,2 4 

--Z +z -4z--Kz + K\n- \-- 

3 2 l-z 4+K 



1 , l-z 
-Kin 

2 z 



+ ln 



K + 4z2 



^ /2z\ ^ /2(l-z) 

— — vKarctan —= + yKarctan z= — 

K + 4(l-z)2 VVk/ V 



K K 



:^z-z +z+-z- -In 



4 l-z 



(100) 



(101) 



B Insertion of emitted particles 

In this appendix we list the relevant formulae for the insertion of one additional four-vector into 
a set of n four-vectors. This insertion satisfies momentum conservation and can be considered 
as the inverse of the (n -I- 1) — n phase space mapping of Catani and Seymour. These insertion 
mappings are also useful for an efficient phase-space integration of the real emission contribution 
in NLO calculations. Therefore we quote in addition the relevant phase space weights. For 
the shower algorithm, these weights are not needed, as they are taken into account through the 
generation of the shower. 



B.l Insertion for final-state particles 

The massless case 

We start with the simplest case, where both the emitter and the spectator are in the final state and 
all particles involved in the dipole splitting are massless. The insertion procedure is identical to 
the one used in [82]. Given the four- vectors pij and together with the three variables y, z and 
we would like to construct pi, pj and pk, such that 

Pi + Pj+Pk = Pij + Pk, pf = pj = Pk = ^- (102) 
In four dimensions we have for the phase space measure 

1 1 271 

dKnres = ^ f (1 -j) f dz J (103) 
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where Sij]^ = {pij + p]^Y = (pi + pj + pk)^. It is convenient to work in the rest frame of P = 
Pij + Pk = Pi + Pj + Pk- We shall orient the frame in such a way, that the spatial components of 
Pk are along the z-direction. When used as a phase space generator we set 



3;= Ml, Z = U2 (|)i = 27lM3, 

where ui,U2 and M3 are three uniformly distributed random numbers in [0, 1]. From 



y = 



Sij + Sik + Sjk' 



z = 



Sik 



Sik + Sjk 



we obtain 



Sij = yP\ Sik = zil-y)P\ sjk = (l-z)(l-y)p2 
If Sij < Sjk we want to have — > pk as Sij — > 0. Define 



Sij + Sjk 



Sij + Sjk 



/Sijk 



Ek^ 



Sik + Sjk 
2y/Sijk ' 



(104) 



(105) 



(106) 



(107) 



Qik = arccos ( 1 



Sik 



2EiEk 



arccos 1 — 



^jk 
lEjEk 



In our coordinate system we have 



p'l = £,(l,sin0,fecos((|),. + ;r),sin0,^sin((|)i, + 7r),cos0,i), 
p'j = Ej(l,smBjkCos^s,smQjkSm(i}s,cosQjk), 



Pk 



£,(1,0,0,1: 



(108) 



(109) 



The momenta p'j, p'- and are related to the momenta pi, pj and pk by a sequence of Lorentz 
transformations back to the original frame 



Pi = ^boosAxy{<^)^xz{^)p'i 



(110) 



and analogously for the other two momenta. The explicit formulae for the Lorentz transforma- 
tions are obtained as follows : Let |P| = a/ (pij + PkY and denote by pk the coordinates of the 
hard momentum pk in the centre of mass system of pij + pk- Pk is given by 



Pk = 



Ep 



Pk-P 



Pk + 



Pk-P 



The angles are then given by 



= arccos 



2EkE'^-2pk-p'^ 



Pk 



\p'k\ 



\P\{Ep + \P\) \P\ 



(|) = arctan ( ^ 

\Pk 



(111) 



(112) 
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For the case considered here particle k is massless and the formula for reduces to 

= arccos 



(113) 



The explicit form of the rotations is 

A.z(e) 



f I \ 

COS0 sin0 

10 

\ -sin0 COS0 / 



The boost p = KhoostQ is given by 



qP 



\P\{Ep+\P\) \P 



/ 1 \ 

cos(t) — sincj) 

sin(|) cos{|) 

\ 1 y 



+ ^\p 



(114) 



(115) 



The weight is given by 



w 



167t2 



(116) 



The massive case 

We now consider the case of final state particles with arbitrary masses: 

~2 22 22 2 -2 2 2 

Pij = ^ip Pi=mi, Pj = mj, p,^ = p^ = m,^. 
The dipole phase space reads [63] 



(117) 



327i3 



y+ z+ 2% 

{l-/4-//j-4f[l{l,/4j,/4)Y''ldy {l-y) JdzJ (118) 



z- 



where 



Sijk= {Pij + Pkf, Ml 



m 



X{x,y,z) = x^ + y^ + z^-2xy-2yz-2zy. (119) 



The integration boundaries are given by 



1 2a'/c(1-A'/J 

y+ = ^--. — .2 ..2 ..2 ' y 



(120) 



d+M] + (}-d-M)-Mk)y 



(121) 
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The general formula for the relative velocities is Vp^q = ^/l — p^q^/{pq)- In our case the relative 
velocities are given by 



Vij,k = 



1 -^/-A/j-A/^J 3-2-4^/2/^2 



(122) 



For the phase space generation we set 

y= {y+-y-)ui+y-, z= {z+ - Z-)u2+Z-, {^j = 27rM3. (123) 

We again work in the rest frame of P = pij + Pk = Pi + Pj + Ph such that the spatial components 
of Pk are along the z-direction: 

Pij={Eij,OA-\Pk\), Pk={Ek,OA\h\)- (124) 
For the invariants we have 

^PiPj = yiP^-mj-mj-ml), 
2piPk = z{\-y){P^-m}-m)-ml), 

2pjPk = {\-z){\-y){P^-m}-m]-ml). (125) 
The invariants are related to y and z as follows: 

'^PiPj _ 2piPk 



y = 



2pipj + 2pipk + 2pjpk' 



z = 



2pipk + 2pjpk 



In our chosen frame 



The energies are obtained from the invariants as follows: 

Ei -- 



Ej = 



Ek 





-2pjpk + mj -mj 










^ijk 


-2piPk-mf + mj 










^ijk 


-2pipj-mj-mj 




2^/^ 



(126) 



p'i = 


Pi KJ^ 
\Pi\ 


sin e,fc cos ((t)i + 7c) , sin 0,^ sm{<ifs + 7c) , cos 6,^) , 




Pi = 




-, sin0y^;cos(^^, sine^-;tsin(^^,cos0;;t), 




p'k = 


\ , Ek 

\Pk\ 


,0,0,1). 


(127) 



(128) 
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For the angles we have 

e,, = arccos ( , 6,, = arccos ( ^^^t'^^fj^' ) ■ (129) 

V ^\Pi\\Pk\ ) \ '^\Pi\\Pk\ ) 

The momenta p\, p'j and are related to the momenta pi, pj and p^ by the same sequence of 
Lorentz transformations as in eq. (110). The weight is 

= {^i^-i4-Mj-Mlf[^ii,i4jA)VHi-y)iy+-y-)iz+-z-). (iso) 

B.2 Insertion for an antenna between an initial-state and a final state 

The massless case 

Here the {n + 1) -particle phase space is given by a convolution: 

1 

d^n+l = j dxd^nixpa) d^dipole- (131) 

The dipole phase space reads: 

d<s?dipoie = ^ 22n^"^ J ^^J ^^^^'^ 



The angle (j), parametrises the solid angle perpendicular to pij and xpa. Therefore we can treat 
the case of a final-state emitter with an initial-state spectator as well as the case of an initial- 
state emitter with a final-state spectator at the same time, x and z are related to the invariants as 
follows: 

^ ^ -^PiPa - 2p jPa - 2pip j „ ^ -2piPa ^^^^^ 

-2piPa - 2p jPa ' -'2-PiPa " 2p jPa 

For the phase space generation we set 

X—l—Ui^ Z — U2, {|)^ = 27lM3. (134) 

We denote Q = pij +xpa = Pi + Pj + Pa- It is convenient to work in the rest frame of P = 
Pi + Pj = Q~ Pa and to orient the frame such that pa is along the z-axis. For the invariants we 
have 

2piPj={-Q^)^, 2piPa = ^Q\ 2pjPa = ^Q\ (135) 



I I 1 271 

\2pijPa\ 
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In this frame 



p'i = £',(l,sin0;acos(|)j,sin0;asin(|)^,cose,a). 
p'j = £,(l,-sine,flCOS(|)^,-sine,aSin{^^,-cos0,, 
Pa = {-\Ea\,0,0,\Ea\sign{pl')). 



la) 7 



We have 



sign(/7^')(-l + 



'^PiPa 



(136) 



(137) 



The momenta p'-, p'j are again related to the momenta pu pj by a sequence of Lorentz transfor- 
mations as in eq. (1 10). The weight is given by 



w 

The massive case 

The dipole phase space now reads: 

d^dipole = 

The integration boundaries are given by 

z+ = 1, 

where 



(138) 



\^PijPa\ 
327t3 
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j dz j d^s 



(139) 



ni: 



l-JC+^2 



xni; 



(140) 



(141) 



\^PijPa\ \Q^-mj\ 

We consider only the case where niQ = mi = m and all other masses are zero. For the phase space 
generation we set 

X=l-Ml, Z= {Z+—Z-)U2+Z-, (l)i = 27tM3. (142) 

For the invariants we have now 

2piPj = (-2' + mj)^, Ipipa = - (e' - m}) , Ipjpa = (e" - m}) . (143) 



We parametrise the momenta as 



Pi 



I A I ( 7^ , sin 0,fl cos , sin 0,a sin , cos ^ ) , 

\Pi\ 



Pj = 1^1(1, -sin0,acos(^,,- sin 0,asin(l)^,-cos0;a), 
Pa = {-\Ea\, 0,0, \Ea\sign{p^J)). 



(144) 
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Then 



Ea = -r^JP- Pa) , Qia = arCCOS 



^,^{2EiEa-lpiPa) 



'^\Pi\{-Ea) 



(145) 



2\P\ ' 

The momenta p'^, p'j are again related to the momenta pu pj by a sequence of Lorentz transfor- 
mations as in eq. (110). The weight is given by 



w = 



167c2jc 

B.3 Insertion for an initial-state antenna 



{Z+-Z-). 



(146) 



Here we only have to consider the case where all particles are massless. In this case we transform 
all the final state momenta. The {n + 1) -particle phase space is given by a convolution: 



1 

d(^n+l = j dx d^n{xpa) d^dipole- 



The dipole phase space reads: 



d^dipole 



l^PaPbl 
327l3 



l-x 2n 



dv J d(^s- 







The variable v is given by 



V = 



-2paPi 
2paPb ' 



For the phase space generation we set 

x=l — ui, V = (1 — x)(l — M2), (|)i = 271113. 



We denote 



We have 



K = -pa -Pb- Pu K = -pai - Pb- 



1 



Pa = -Pai, 
X 

Pi = A^ooif^/(l,sin0,acos(^^,sin0,asin{^^,cos0,a) 
Pb = Pb, 

with Ei and Qia given in the rest frame of pa + Pb by 



17 1 A r K^-'^PaPb Q 

= - 7: V^PaPb , Ei = — , eia = arccos 



sign(Pa) -1 + 



2piPa 
2EiEa 



(147) 



(148) 



(149) 



(150) 



(151) 



(152) 



(153) 
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pa denotes pa in the rest frame of pa + Pb- ^boost transforms from the rest frame of pa + Pb to 
the lab frame. All other final state momenta are transformed with 

A-i = gMV_2\ ^ \ ) "rl—^. (154) 

The weight is given by 
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Figure 1: The first moments of the thrust distribution, the C-parameter distribution and the D- 
parameter distribution. The results are from the parton shower for two different values of the 

cut-off scale Qmin- 
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Figure 2: The distributions for the four-jet angles. From top-left to bottom-right: The modified 
Nachtmann-Reiter angle, the Komer-Schierholz-Willrodt angle, the Bengtsson-Zerwas angle and 
the angle a34 between the smallest energy jets. As cut-off parameter Q,nin = 1 GeV is used. 
Shown are the result from the "strict leading-colour approximation" and the "modified leading- 
colour approximation". 
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Figure 3: The transverse momentum distribution and the rapidity distribution of the lepton pair 
for Z/t* -production for the Tevatron and the LHC. As cut-off parameter Q^m = 1 GeV is used. 
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